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\/y ■ Abstract 

' Considering the coupled envelope equations in nonlinear couplers, the question of in- 

tegrability is attempted. It is explicitly shown that Hirota's bilinear method is one of 
the simple and alternative techniques to Painleve analysis to obtain the integrability 
conditions of the coupled nonlinear Schrodinger (CNLS) type equations. We also show 
that the coupled Hirota equation introduced by Tasgal and Potasek is the next hier- 
archy of the inverse scattering solvable CNLS equation. The results are in agreement 

(N . with the known results. 

O 

^ ■ 1 Introduction 

00 

The propagation of temporal soliton envelopes in nonlinear optical media has been pre- 
dicted and demonstrated experimentally [Ij— 1|. This prediction arises from the opportu- 
nity of reducing the Maxwell equations which govern the propagation to a single, com- 
pletely integrable, partial differential equation(PDE) where one is left with the time and 
the propagation distance as independent variables [Q] . The best and well studied example 
^ ■ is the nonlinear Schrodinger (NLS) equation 

iq z + q tt + 2 | q \ 2 q = (1) 

which describes the wave propagation of picosecond-pulse envelopes q(z, t) in a lossless 
single-mode fibre. Eq.(l) is one of the completely integrable iV-soliton possessing systems. 
Eventhough eq.(l) may adequately describe the propagation in a single-mode waveguides, 
routing-, switching operations, and other optical effects which involve soliton pulses re- 
quire the interaction between two or more modes Q]— §• In general, the coupled mode 
approach M still permits description of the pulse propagation in a multi-mode waveguide 
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by means of vector versions of eq.(l). Although these systems of equations are no longer 
integrable, one may obtain quantitative information about the pulse propagation by re- 
sorting to numerical or perturbative methods. Situations of physical interest that can be 
described by CNLS equations include two parallel waveguides coupled through evansecnt 
field overlap, the coupling of two polarizations modes in uniform guides, and so on. Note 
that the study of the propagation of optical solitons in multi-mode nonlinear couplers, 
besides being important from the theoretical point of view, is also important in view of 
their possible applications. Recent advances have permitted proof that solitons are ideal 
candidates for performing all-optical switching operations in nonlinear couplers. In fact, 
their stability leads to the possibility of controlling the coupling of the whole pulse, by 
means of changing the input power of a single pulse, or the phase of the superimposed 
pulse. 

For the past few years, the wave propagation in nonlinear couplers have been inves- 
tigated from theoretical, numerical and experimental point of view 

@|-fnj. From the 

detailed theoretical investigations, several integrable soliton possessing coupled NLS equa- 
tions have been introduced || [9|, 1C, 12, 16, IE]. The experimental results reveals that 



these integrable soliton models are not accurate enough to explain the wave propagation 
in nonlinear couplers. On the otherhand, theoretically several new results and concepts 
have been introduced. Thus it is clear that to propagate exact soliton through nonlinear 
couplers, one has to look for different conditions among the parameters involved in the 
system. In this letter, we present new type of bilinear conditions, which may be very 
useful for the experimental realization of solitons in nonlinear couplers. 

It is well known that Painleve singularity structure analysis is one of the systematic and 
powerful method in nonlinear science to identify the integrability conditions of nonlinear 
partial differential equations (NPDEs) In recent years, this method has been ap- 

plied to a very large number of NPDEs and also systematically estabilished the complete 
integrability properties like Lax pair, Backhand transformation, bilinear transformation, 
soliton solutions, and so on [13, 20, [Tj], 22]. For CNLS type equations, this method is found 
to be cumbersome and the construction of the Lax pair from this analysis is still an open 
question. Also, after getting the necessary condition(s) for integrability, one has to apply 
bilinear or other methods to prove the sufficient condition for the complete integrability. 
The main aim of this letter is to show explicitly that Hirota's bilinear approach [23] can 
be used as an alternative method to Painleve analysis to identify the complete and partial 
integrability conditions of CNLS equations with many parameters. Using this method, we 
show that one can obtain the more general bilinear (integrability) conditions of coupled 
nonlinear equations in a very simple way. 



2 Two Coupled Nonlinear Schrodinger Equation 

There are several methods to derive a set of CNLS equations, depending on the physical 
situation that is being modeled. A fairly general and frequently studied CNLS equation 
is of the form |, § §, §, 0, 0, 0] 

iu z = c\u t t + 2 (a | u | 2 +j3 \ v | 2 ) u, (2) 
iv z = c 2 v a + 2 {[3 | u | 2 +7 | v | 2 ) v. (3) 
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Eqs.(2-3) are shown to be completely integarble for the following two cases: 

(0 : ci = c 2 , a = (3 = 7; (ii) : c\ = -c 2 , a = -(5 = 7. (4) 

To obtain the (bilinear) integr ability conditions and soliton solutions, we introduce the 
following transformation |2^] 

u=-, v = -, (5) 

where G and H are complex functions and F is a real function and the bilinear operator 
is defined by 

D " D " = (I " ikT {I ~ £)"«<*.'W) 1 - w 

Substituting eq.(5) in eqs(2-3), we obtain 

(iD z - Cl D 2 t ) G-F = 0, (iD z - c 2 D 2 ) H ■ F = 0, (7) 

D 2 t F ■ F = — (aGG* + 0HH*) , D\F -F = — (jHH* + f3GG* ) . (8) 

Cl C2 

So the left hand sides of eq.(8) become equal. Hence the right hand sides of eq.(8) should 
also be equal which is true only when 

c 2 a = ci(3, c 2 P = C17. (9) 

The above conditions can be obtained by equating the coefficients of GG* and HH* 
respectively in eq.(8). One can easily check that eq.(9) admits the already known integra- 
bility conditions (4). Solving eq.(9), we obtain 

c 2 = ± A /^ci, ai = (5 2 . (10) 
V a 

It should be noted that conditions (10) are more general than (4). As the bilinear condi- 
tions (9) or (10) are in terms of all parameters, we believe that these conditions may be 
very useful for the experimental generation of solitons in nonlinear couplers. 

Having obtained the bilinear conditions and bilinear forms, our next aim is to obtain 
the soliton solutions. Now to generate the soliton solutions, we expand the dependent 
variables in terms of e and the solutions can be constructed as usual. For instance, if 
G = sGi, H = eH\ and F = 1 + e 2 F 2 , the solutions are found to be 

Gl = e m H = e*+C, F 2 = El + ^ , (11) 

ci(«i + k* x y 

where 771 = k\(t — ic\k\z) + rj 2 = ki(t — ic 2 k\z) + rf( and £ is a complex constant. From 
the detailed investigations, we find that eqs.(2-3) admit the higher order soliton solutions 
for the conditions given in eq.(4) and the resultant solutions are in agreement with the 
resluts reported earlier (9|, |TJ, |l~5|| . So like fixing the parameters in the arbitrary analysis 
of Painleve method, using bilinear method one can also get the complete integrability 
conditions from the generation of higher order soliton solutions. 
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3 Three Coupled NLS Equations 

Considering the systems of three coupled NLS equations of the form [14] 

iu z = c\Uu + 2 (a | u \ 2 +3 \ v | 2 +5 j w | 2 ) u, (12) 



iv z = c 2 v tt + 2 (3 | u | 2 +7 | v | 2 +r | w | 2 ) v, (13) 

= c 3 w u + 2(5 \ u\ 2 +T \ v \ 2 +A | u; | 2 ) w. (14) 

Using the similar procedure the bilinear conditions are found to be of the form: 

c 2 c 3 a = cic 3 (3 = cic 2 5, c 2 c 3 8 = C1C37 = cic 2 T, c 2 c 3 5 = C1C37 = cic 2 A. (15) 

Here also we find that eqs.(15) satisfy the integrability conditions reported through 
Painleve analysis jl4| and the soliton solutions can also be constructed. 



4 Coupled Hirota equation 



In this section, we will discuss the integrability of the coupled Hirota equation [16]. In 
order to increase the bit rates it is necessary to decrease the pulse width of the order of 
femtosecond. As pulse lengths become comparable to the wavelength, however, eqs.(2-3) 
are inadequate, as additional effects must now be considered. By considering the above 



facts, a generalized coupled Hirota equation reads as |16] 
iu z + c\Uu + 2 (a I u | 2 +3 I v I 2 ) u 

(16) 

-ie [u m + (2/ii I u | 2 +v\ I v | 2 ) u t + v 1 uv*v t ] = 0, 
iv z + c 2 v u + 2 (3 I u | 2 +7 I v | 2 ) v 

(17) 

-ie [vut + \V2 I u I 2 +2^2 I v \ 2 ) v t + v 2 u*vut\ = 0. 

When ci = c 2 = — , a = 3 = 7 = A, [i\ = v\ = [i 2 = v 2 = 3, eqs. (16-17) take the 

form of the equations in [O]. In eqs. (16-17), we have not included the linear effects, 
which can always be removed through suitable transformation in u and v [15| and, for our 



convenience, we introduce the new parameters v 2 and \x 2 . The bilinear forms of (16-17) 
are in the form 

(iD z + ciD 2 - ieDf) G ■ F = 0, (iD z + c 2 D 2 - ieD^) H ■ F = 0, (18) 

D 2 F ■ F = — (aGG* + BHH*) , D 2 F ■ F = - (mGG* + ^HH*) , (19) 

ci 3 

D\F ■ F = — (jHH* + BGG*) , D 2 F ■ F = -(u 2 GG* + fi 2 HH*) , (20) 

c 2 3 



The more general bilinear conditions are: 

ci 3 c 2 3 ' ci 3 c 2 3 



(21) 
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Here also we find that eq.(21) admits the higher order soliton solutions for the conditions 



ci = C2 , a = (3 = 7, and [i\ = v\ = ^2 = ^2 = 3 jig, 17]. Exact iV-envelope soliton 
solutions of eqs. (16-17) can also be expressed in a similar way. It is interesting to note 
that, like CNLS, eqs. (16-17) also found to be completely integrable for another conditions 
ci = — C2 , a = — /3 = 7, and fi± = —v\ = = ^2 = 3. It should be noted that 
the latter condition is one of the new integrable system obtained through this method. 
Recently, using multi-Hamiltonian formalism, we also verified the above complete integra- 
bility conditions and found that the above integrable cases are the next hierarchy of the 



CNLS equations with conditions (4) [18]. Hence, the coupled Hirota equation introduced 



by Tasgal and Potasek [16[] is the next hierarchy of the 1ST solvable CNLS equation-I. 



5 Coupled Higher Order NLS Equation 



We consider the coupled higher order NLS equation of the form [19] 
iu z + c\Uu + 2 (a | u | 2 +/3 \ v | 2 ) u 



(22) 

is [u ttt + 6(\ u | 2 + | v | 2 ) u t + 3(| u | 2 + | v | 2 Lu| 



>t 



iv z + c 2 v t t + 2 {(3 | u | 2 +7 | v | 2 ) v 

(23) 

-is [v ttt + 6(\ u | 2 + | v | 2 ) v t + 3(| u | 2 + | v \ 2 ) t v] = 0. 

Equations (22-23) can be derived from the higher order NLS ||] by considering the 
electromagnetic wave £ as a sum of left and right polarized waves |19|| . When c\ = C2 = — 
and a = f3 = 7 = 1, the bilinear forms are found to be 



iD z + -D 2 t - ieB$ ) G ■ F = 0, ( iD z + -D 2 - ieD* ) H ■ F = 0, (24) 



±D 2 -ieD^G-F = 0, (iD, + ± 

D 2 F ■ F = 4(GG* + HH*), (25) 
D t G-G* = 0, D t G-H* = 0, D t G* ■ H = 0, D t G ■ H = 0, D t H ■ H* = 0.(26) 
The one soliton solution of eqs. (22-23) is constructed as 

u = aikl sech (kit + ek\z + ry?) exp ( % _hll\ (27) 

V2(| ai | 2 + I a 2 | 2 ) 1 1 ' lj \ 2 j' 1 > 

v = a2kl sech (kit + ek\z + 77?) exp ( l Jlll\ (28) 

yffl ai | 2 + I a 2 | 2 ) 11 1 /U V 2 ; y ' 

where a\ and a 2 are integration constants. 

Thus in this letter, we have shown that Hirota's bilinear method is one of the alternative 
formalism to Painleve singularity structure analysis to identify the integrability of coupled 
nonlinear Schrodinger equations with many parameters. This analysis is found to be 
very simple when compared with the Painleve analysis. Also from the optical point of 
view, the general bilinear conditions obtained through bilinear form may be useful for 
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the experimental generation of solitons in nonlinear couplers. Another advantage of this 
method is that one can, in a simple manner, construct the soliton solutions for all integrable 
cases. We have also pointed out that the coupled Hirota equation studied by Tasgal 
and Potasek is the next hierarchy of the 1ST solvable CNLS equation. We found that 
the conditions obtained through this method are in agreement with the results reported 
earlier. It will be interesting to investigate the nature of the solutions (at least two soliton 
solutions) for the new integrable cases. Work is in progress in this direction. 
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